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ABSTRACT 

We study the Extended Chandra Deep Field South (E-CDFS) Very Large Array 
sample, which reaches a flux density limit at 1.4 GHz of 32.5 pdy at the field 
centre and redshift ~ 4, and covers ~ 0.3 deg 2 . Number counts are presented for 
the whole sample while the evolutionary properties and luminosity functions are 
derived for active galactic nuclei (AGN). The faint radio sky contains two totally 
distinct AGN populations, characterised by very different evolutions, luminosity 
functions, and Eddington ratios: radio-quiet (RQ)/radiative-mode, and radio- 
loud/jet-mode AGN. The radio power of RQ AGN evolves oc (1 + z) 2 ' 5 , similarly 
to star-forming galaxies, while the number density of radio-loud ones has a peak 
at z ~ 0.5 and then declines at higher redshifts. The number density of radio- 
selected RQ AGN is consistent with that of X-ray selected AGN, which shows 
that we are sampling the same population. The unbiased fraction of radiative¬ 
mode RL AGN, derived from our own and previously published data, is a strong 
function of radio power, decreasing from ~ 0.5 at Pi. 4 GHz ~ 10 24 W Hz -1 to 
~ 0.04 at Pi. 4 GHz ~ 10 22 W Hz -1 . Thanks to our enlarged sample, which now 
includes ~ 700 radio sources, we also confirm and strengthen our previous results 
on the source population of the faint radio sky: star-forming galaxies start to 
dominate the radio sky only below ~ 0.1 mJy, which is also where radio-quiet 
AGN overtake radio-loud ones. 

Key words: galaxies: active — galaxies: evolution — galaxies: starburst 
radio continuum: galaxies — quasars: general - surveys 


1 INTRODUCTION 

Soon after the discovery of quasars (Schmidt 1963) it was 
realized that the majority of them were not as strong ra¬ 
dio sources as the first quasars, as they were undetected 
by the radio telescopes of the time: they were “radio¬ 
quiet” (Sandage 1965). It was later understood that these 
sources were actually only “radio-faint”, as for the same 
optical power their radio powers were ss 3 — 4 orders 
of magnitude smaller than their radio-loud (RL) coun¬ 
terparts, but the name stuck. Radio-quiet (RQ) active 
galactic nuclei (AGN) were until recently normally found 
in optically selected samples and are characterized by rel¬ 
atively low radio-to-optical flux density ratios (R < 10) 
and radio powers (Pi.4GHz < 10 24 W Hz -1 locally). 

Innumerable studies have compared the properties 
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of the two AGN classes in various bands to try and shed 
light on their inherent differences. As a result, the dis¬ 
tinction between the two types of AGN has turned out 
to be not simply a matter of semantics: the two classes 
represent intrinsically different objects, with RL AGN 
emitting a large fraction of their energy non-thermally 
and in association with powerful relativistic jets, while 
the multi-wavelength emission of RQ AGN is dominated 
by thermal emission, directly or indirectly related to the 
accretion disk. The host galaxies are also different, with 
those of RL AGN being elliptical while those of RQ ones, 
excluding the most powerful ones, being spiral (e.g. Dun¬ 
lop et al. 2003). 

However, fifty years after the discovery of quasars the 
question “Why do only a minority of galaxies that contain 
an AGN have jets?” is still unanswered. One problem is 
that, while it would be important to compare the proper¬ 
ties of the two AGN classes in the band where they differ 
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most (i.e. the radio band) this requires the identification 
of AGN in deep (< 1 mjy) radio fields, as RQ AGN are, 
by definition, radio faint. This has been possible only re¬ 
cently (e.g. Padovani et al. 2009, 2011; Simpson et al. 
2012; White et al. 2015, and references therein). 

Our group started addressing this topic, together 
with others including the broader issue of the source pop¬ 
ulation in very deep radio fields, initially in the Chandra 
Deep Field South (CDFS). This was done by defining a 
complete sample of 198 radio sources sources reaching 
~ 43 /rJy over 0.2 deg 2 at 1.4 GHz using the National 
Radio Astronomy Observatory (NRAO) Very Large Ar¬ 
ray (VLA) through a multi-wavelength approach (Keller- 
mann et al. 2008; Mainieri et al. 2008; Tozzi et al. 2009; 
Padovani et al. 2009, 2011). 

Miller et al. (2008, 2013) expanded on this by ob¬ 
serving the so-called Extended CDFS (E-CDFS), again 
using the VLA, down to ~ 30 /iJy at 5a, in a 2.8” x 1.6” 
beam over ~ 0.3 deg 2 . This resulted in a sample of almost 
900 sources. We have started exploiting these new radio 
data with the aim of addressing the issues of the faint 
radio source population and RQ AGN in more detail, 
given the larger and slightly deeper E-CDFS sample, as 
compared to the CDFS one. In particular, Bonzini et al. 
(2012) have identified the optical and infrared (IR) coun¬ 
terparts of the E-CDFS sources, finding reliable matches 
and redshifts for ~ 95% and ~ 81% of them respectively, 
while Vattakunnel et al. (2012) have identified the X-ray 
counterparts and studied the radio - X-ray correlation for 
star forming galaxies (SFG). Finally, Bonzini et al. (2013) 
have provided reliable source classification and studied 
the host galaxy properties of RQ AGN, RL AGN, and 
SFG. 

The main aims of this paper are: 

(i) present the most accurate number counts for the 
various classes of sub-mjy sources down to ~ 30 /rjy; 

(ii) study the evolution and luminosity functions (LFs) 
of sub-mjy AGN. In particular, our CDFS paper on this 
topic (Padovani et al. 2011) was affected by small num¬ 
ber statistics, particularly for RQ AGN, and by the large 
fraction (~ 50%) of upper limits on the far-IR-to-radio 
flux density ratio ( q ) used for classifying sources. As dis¬ 
cussed at length in Padovani et al. (2011), this topic has 
important implications for constraining the mechanism 
behind radio emission in RQ AGN and allow a proper 
comparison with RL AGN; 

(iii) investigate more deeply the density evolution of 
RL AGN found in Padovani et al. (2011) (oc (1 + 2 )” 1 ' 8 ); 

(iv) provide a new approach to the issue of the fraction 
of RL sources within the AGN population. 

Items (i) and (ii) are very relevant, together with the 
study of the evolution of SFG, also to the issue of the link 
between the black holes at the centre of AGN and their 
host galaxies. Moreover, they are also critical to predict 
the source population at radio flux densities < 1 /rjy, 
which are important, for example, for the Square Kilo¬ 
metre Array (SKA). All existing estimates, in fact, had to 
rely, for obvious reasons, on extrapolations and are based 
on high flux density samples. This affects particularly the 
highest redshifts, which can better be probed at fainter 
flux densities. 


Our physical definition of the various classes of sub- 
mjy sources follows that of Padovani et al. (2011), to 
which we refer the reader. Following Heckman & Best 
(2014) we will also use the terms “radiative-mode” and 
“jet-mode” AGN 1 . In short, radiative-mode RL AGN 
include radio quasars and high-excitation radio galax¬ 
ies, while jet-mode RL AGN refer to low-excitation ra¬ 
dio galaxies. There is considerable overlap between jet¬ 
mode AGN and the radio sources morphologically clas¬ 
sified as Fanaroff-Riley type I (FR I) (Fanaroff & Riley 
1974) and also between RL radiative-mode AGN and FR 
IIs, although there is a sizeable population of jet-mode 
FR IIs and a smaller one of radiative-mode FR I (e.g. 
Gendre et al. 2013, and references therein). As for RQ 
AGN, radiative-mode sources are the “classical” broad- 
and narrow-lined AGN (type I and II), while the jet¬ 
mode ones are the so-called liners. Note that the two 
classes have also widely different Eddington ratios, with 
radiative-mode and jet-mode sources respectively above 
and below L/L Edd ~ 0.01. 

In this paper we present the general properties (num¬ 
ber counts, redshift distribution, etc.) of the whole sam¬ 
ple and then concentrate on the evolution and LFs of 
AGN. We will discuss in detail SFG in a future publica¬ 
tion (Padovani et al., in preparation). 

Section 2 describes the VLA-E-CDFS sample, while 
Section 3 presents the so-far deepest determination of 
sub-mjy number counts by class. Section 4 studies 
the sample evolution while Section 5 derives the AGN 
LFs and their evolution. Finally, Section 6 discusses 
our results, while section 7 summarizes our conclusions. 
Throughout this paper spectral indices are defined by 
S„ oc v~ a and the values Ho = 70 km s -1 Mpc -1 , 
Om = 0.3, and IIa = 0.7 have been used. 

2 THE SAMPLE 
2.1 Redshifts 

Our original sample is described in Bonzini et al. (2012) 
and Miller et al. (2013) and includes 883 radio sources 
detected in 0.324 deg 2 at 1.4 GHz in a deep VLA sur¬ 
vey of the E-CDFS. The average 5a flux density limit 
is ~ 37 fiJy, reaching ~ 30 /rJy in the field centre. The 
fraction of sources with redshift information is ~ 81%, 
which is not ideal for determining the evolutionary prop¬ 
erties and LFs of our sample. We have then excluded 
the outermost part of the field, where we do not have 
enough ancillary multi-wavelength data to provide rea¬ 
sonable photometric redshifts. The resulting sub-sample 
covers an area of 0.285 deg 2 (see Bonzini et al. 2013, for 
details) and includes 765 radio sources with radio flux 
density ^ 32.5 /rJy, 87% of which have redshifts (of these, 
~ 40% are spectroscopic and ~ 60% are photometric). 
This is the sample we use to build the number counts 
in Sect. 3. To further increase the fraction of sources 
with redshift, we have made another cut and considered 
only sources with a 3.6/rm flux density / 3 . 6 /jm > 1 mi- 
crojy, which is approximately the completeness limit of 

1 See Table 4 of Heckman & Best (2014). 
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the IRAC observations at 3.6/rm (only four sources have 
/ 3 . 6 , 1 m < 1 microjy). Coupled with a radio flux density 
limit of 32.5 /rjy, this provides us with clean selection 
criteria and with a sample of 680 radio sources, 92% of 
which have redshift (the same fraction as in Padovani 
et al. 2011). We estimate the redshift for the ~ 8 % of 
the objects in the sample without observed redshift from 
their / 3 . 6 , 1 m, as detailed in Appendix A. To further min¬ 
imise the effects of these redshift estimates we also define 
two subsamples: AGN (and RQ AGN) with z 4 3.66 and 
SFG with z < 3.25. Above these two values, in fact, 80% 
and 95% of the redshifts, respectively, are estimated. By 
imposing such cuts, we are then only assuming that the 
excluded sources are at high redshifts, which makes sense 
based on their / 3 . 6 , 1 m values (see Appendix A), without 
actually using their values. This is the sample we use 
mostly to study the evolution and LF of AGN in Sect. 4 
and 5. 

Note that the effect of these redshift estimates on our 
results is minimal, for two reasons: 1. our final redshift 
incompleteness is very small (~ 3%, ~ 5%, and ~ 1.5% 
for all AGN, RQ, and RL AGN respectively: see Sect. 
4.1); 2. redshift affects 14/14 values (Sect. 4.1) much less 
than flux density. Nevertheless, error bars on the LFs 
are evaluated using the number of sources per bin with 
redshift determination only and when binning in redshift 
the percentages of estimated redshifts are also given for 
each bin. 


2.2 Classification 

Source classification has been done by Bonzini et al. 
(2013), expanding upon the scheme of Padovani et al. 
(2011). In short, by using 52406 s, that is the logarithm 
of the ratio between the observed 24 /im and 1.4 GHz 
flux densities, we could define an “SFGs locus” based 
on the radio - FIR. correlation for SFG. Sources below 
this locus display a radio excess that is the signature of 
an AGN contribution to the radio luminosity and were 
therefore classified as RL AGN. Within and above this 
locus, a source was classified as a RQ AGN if there was 
evidence of AGN activity in the other bands considered. 
Namely, if it had a hard band (2 — 10 keV) X-ray luminos¬ 
ity greater than 10 42 erg s _1 or it lay within the “AGN 
wedge” of the IRAC colour-colour diagram, as defined by 
Donley et al. (2012). Otherwise, the object was classified 
as an SFG. We point out that our AGN classification 
includes all sources independently of optical appearance 
(i.e., quasar-like and galaxy-like). We refer the reader to 
Bonzini et al. (2013) for further details. 

The issue of the classification of our sources into 
radiative- or jet-mode sources has been briefly discussed 
by Bonzini et al. (2013). The optimal approach to be able 
to make such a distinction requires rest-frame optical/UV 
emission line ratios but, given the relatively high redshift 
of our sample and the incomplete spectroscopic coverage, 
this is only possible for a handful of objects. Bonzini et 
al. (2013), using a criterion based on the 22 /im power 
proposed by Giirkan, Hardcastle, & Jarvis (2014), esti¬ 
mated that the large majority of our RL sources are of 
the jet-mode type; this is confirmed by their LF (Sect. 


5.2.1). As regards our RQ AGN, given our selection cri¬ 
teria, which include also a cut in X-ray power, we expect 
them to be for the large part of the radiative-mode type. 
This issue is addressed in Sect. 6.2. 

We stress that ours is the deepest sample for which 
both the evolution and LFs of radio sources has been 
studied. 


3 THE SUB-MJY NUMBER COUNTS 

More than half a century ago, the earliest radio source 
counts gave the first indications of cosmic evolution (Ryle 
& Scheuer 1955) but the results were not without con¬ 
troversy. The central problem in understanding source 
counts is the conflict between source confusion at lower 
resolution and missing lower surface brightness sources 
or underestimating their flux density at high resolution. 
These issues, while better understood, still exist today. 
There is no “correct” resolution to use that avoids the 
effects of confusion and resolution, but there is an opti¬ 
mum resolution depending on the areal density of sources. 
Modern deep radio surveys reach flux density levels of 
only a few pjy where the corresponding source density 
is ~ 5 sources arcmhU 2 at S = 10 pjy (Vernstrom et 
al. 2015). So a resolution of at least a few arc seconds 
is needed to unambiguously separate sources. But many 
/rjy sources are typically an arcsecond in size, so are par¬ 
tially resolved with a one or two arcsecond beam (e.g. 
Muxlow et al. 2005). 

Published deep radio source counts are widely dis¬ 
crepant, even when made with the same instrument, and 
even with different investigators using the same instru¬ 
ment in the same field (Condon et al. 2012). Since source 
detection depends on the peak flux density, the effects 
of resolution cause two problems. First, if a source is re¬ 
solved the measured peak flux density needs to be multi¬ 
plied by the ratio of integrated to peak flux density. For 
strong sources, where both the measured and peak flux 
densities can be accurately determined this is straight¬ 
forward. But for sources with lower signal-to-noise ratio 
(SNR) Gaussian fitting routines give systematically spu¬ 
rious large sizes and a corresponding flux density bias 
(e.g. Condon 1997). Even more complicated are sources 
where the peak flux density is below the detection limit, 
but with integrated flux densities which are above the 
detection limit. To account for these non-detections one 
normally may use a source size distribution equivalent to 
that of stronger survey sources, but this does not neces¬ 
sarily give the correct answer. As described by Miller et 
al. (2013) we have used peak flux densities, except where 
we measure a significantly higher integrated value, a pro¬ 
cedure which depends on the local SNR. Because of the 
then limitations of the VLA at the time of the Miller et al. 
(2013) survey, the 25 MHz IF bands could be divided into 
only seven sub bands each of 3.125 MHz which resulted in 
significant bandwidth smoothing away from the field cen¬ 
tre. Because each source may appear in more than one 
subfield field, the corrections for bandwidth smoothing 
are not straightforward (Miller et al. 2013). 

An additional problem arises when one physical 
source is split into two or more distinct components. 
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S (1.4 GHz) [mJy] S (1.4 GHz) [mJy] 

Figure 1. a) (left) the Euclidean normalised 1.4 GHz ECDFS source counts (black filled squares) compared to the source counts 
at 1.4 GHz from the compilation of Condon (1984) (red filled circles); b) (right) the Euclidean normalised 1.4 GHz ECDFS source 
counts for the whole sample (black filled squares) and the various classes of radio sources: SFGs (green diamonds), all AGN 
(magenta triangles), RQ AGN (blue circles), and RL AGN (red squares). Error bars correspond to la Poisson errors (Gehrels 
1986). See text for details. 


Table 1. Euclidean normalised 1.4 GHz counts. 
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Unless there is a connection between the radio features, 
there is no way to a-priori distinguish between a rela¬ 
tively strong multi-component source or multiple weaker 
sources. But the derived source count will depend dra¬ 
matically on which interpretation is used. We have ad¬ 
dressed this issue by reference to optical and/or X-ray 
counterparts using a maximum likelihood criteria to de¬ 
cide which optical counterparts are associated with which 
radio features (Bonzini et al. 2012). 

Finally, we note the problem caused by the so-called 
Eddington bias. Because the source counts are steep, that 
is there are more weak sources than strong sources, ran¬ 
dom noise fluctuations cause more weak sources to appear 
above any flux density value, than strong sources which 
fall below this value. As described by Wall & Jenkins 
(2012) the determination of the Eddington bias to the 
observed source count is not straightforward. Following 
the procedure described by Wall & Jenkins (2012), Wall 
(private communication) has calculated the Eddington 


bias corrections to the source count assuming a Euclidean 
slope. We then corrected the number of sources by a fac¬ 
tor, which depended on SNR and ranged from 1/1.194 
(5 < SNR < 6) to 1/1.052 (9 < SNR < 10). These same 
factors were also applied when deriving the LFs. 

To investigate the relative contribution of the differ¬ 
ent source types to the radio population as a function of 
flux density, we consider here the sub-sample within an 
area of 0.285 deg 2 , which provides us with more infor¬ 
mation for this purpose (i.e., we do not make any 3.6/rm 
flux cut). The sensitivity of our survey is a function of 
the position in the field of view, although a much less 
strong one than for the single pointing CDFS sample. In 
more detail, the area covered is equal to 0.285 deg 2 and 
independent of flux density above ~ 60 /rjy, decreases 
quite slowly to 0.279 deg 2 and 0.226 deg 2 at ~ 50 /iJy 
and ~ 40 /rJy respectively, and reaches 0.05 deg 2 at the 
flux density limit (see Miller et al. 2013). We included 
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sources down to a root mean square (rms) noise of 6.5 
pjy, that is a flux density limit of 32.5 pjy. 

Fig. la compares the E-CDFS number count with 
previous determinations at relatively high flux densities 
compiled by Condon (1984). The strong source count 
above ~ 1 Jy shows the familiar steep slope resulting from 
the well know dramatic evolution of radio loud quasars 
and radio galaxies. Below this value, all of the strongest 
sources in Universe have already been included, and the 
count begins to converge. Below 1 mjy, the count steep¬ 
ens (flattens in the normalised plot) again due to the 
emergence of the low luminosity population of RQ AGN 
and SFG, which are discussed further in this paper and 
in the next one of this series. 

Table 1 and Fig. lb 2 present the Euclidean nor¬ 
malised number counts for the various classes. The nor¬ 
malised number counts agree very well in the region of 
overlap (flux densities > 43 pjy) with those of the CDFS 
(Kellermann et al. 2008) 3 and appear to be slightly higher 
at fainter flux densities. The comparison with the num¬ 
ber counts of Padovani et al. (2011) is not shown for clar¬ 
ity but the two determinations agree very well, typically 
within ~ la. The main results of Padovani et al. (2011) 
are also confirmed: AGN dominate at large flux densities 
(> 1 mjy) but SFGs become the dominant population 
below ss 0.1 mjy. Similarly, RL AGN are the predomi¬ 
nant type above 0.1 mjy but their counts drop fast at 
lower values. 

In more detail, AGN make up 43 ± 3% (where the 
errors are based on binomial statistics: Gehrels 1986) of 
sub-mjy sources and their counts are seen to drop at 
lower flux densities, going from 100% of the total at ~ 
10 mjy down to 38% at the survey limit. SFG, on the 
other hand, which represent 57 ± 4% of the sample, are 
missing at high flux densities but become the dominant 
population below ss 0.1 mjy, reaching 62% at the survey 
limit. RQ AGN represent 26±2% (or 61% of all AGN) of 
sub-mjy sources but their fraction appears to increase at 
lower flux densities, where they make up 75% of all AGN 
and « 29% of all sources at the survey limit, up from 
S3 6% respectively at ss 1 mjy. The relative fractions of 
the various classes in our sample are shown in Fig. 6 of 
Bonzini et al. (2013). 

4 EVOLUTION 
4.1 14/14 analysis 

We first study the evolutionary properties of the VLA- 
E-CDFS sample through a variation of the V/V m3bX test 
(Schmidt 1968), the 14/14 test (Avni & Bahcall 1980; 
Morris et al. 1991), that is the ratio between enclosed 
and available volume, since we do not have a single flux 
density limit (see Sect. 3). More details on this and other 


2 A preliminary version of this figure has been shown by 
Padovani et al. (2014). 

3 Due to a small numerical error, the number counts in Kellcr- 
mann et al. (2008) and, consequently, also those in Padovani et 
al. (2009) and Padovani et al. (2011), were too high by ~ 5%. 



redshift 


Figure 2. Fractional redshift distributions for the different 
classes of sources, deconvolved with the appropriate sky cov¬ 
erage. The dashed areas denote redshifts estimated from the 
3.6pm flux density. Error bars represent the 1 a range based 
on Poisson statistics. The dotted lines take into account the 
presence of LSS as discussed in Appendix B. 


statistical techniques used in this paper can be found in 
Padovani et al. (2011). 

We have computed 14/14 values for our sources tak¬ 
ing into account our double flux density limits and the 
appropriate sky coverage (see eqs. 42 and 43 of Avni & 
Bahcall 1980). The IR k-correction was done by using 
the measured near-IR spectral indices. In the radio band 
spectral indices between 1.4 and 4.86 GHz are available 
for ~ 29% of the sources (Kellermann et al. 2008; Huynh 
et al. 2012), while for the rest the mean values for the rel¬ 
evant classes were assumed. Statistical errors are given by 
a = 1/V12 N (Avni & Bahcall 1980). We estimate the 
significance of the deviation from the non-evolutionary 
case by deriving the p -value that the 14/14 distribution 
is similar to a uniform one according to a Kolmogorov- 
Smirnov (KS) test. Similar results are obtained by using 
the deviation from 0.5 (the value expected for no evolu¬ 
tion) of (14/14). To have an initial simple estimate of the 
sample evolution we have also derived the best fit param¬ 
eter fcz, assuming a pure luminosity evolution (PLE) of 
the type P(z) = P(0)(l + z) kL or a pure density evolution 
(PDE) of the type 4>(,z) = 4>(0)(1 + z ) k °, where <f>( 2 ) is 
the luminosity function. 

Our results are shown in Table 2, which gives the 
sample in column (1), the number of sources in column 
(2), the mean redshift in column (3), the percentage of 
sources with redshift estimated from the magnitude in 
column (4), (14/14) in column (5), the p -value that the 
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Table 2. Sample Evolutionary Properties: 14/14, analysis. 


Sample 

N 

<*) 

est. 2 % 

(Ve/Va) 

p- value 

k L a 

k D b 

All sources 

680 

1.35 ± 0.05 

8.4% 

0.609 ±0.011 

< 0.001 

1.6 ±0.1 


Star-forming galaxies 

372 

1.16 ±0.05 

8.9% 

0.656 ±0.015 

< 0.001 

2.1 ± 0.1 


Star-forming galaxies, 2 : ^ 3.25 

356 

1.01 ± 0.04 

5.1% 

0.651 ±0.015 

< 0.001 

2.4 ±0.2 


All AGN 

308 

1.54 ±0.08 

7.8% 

0.550 ±0.016 

< 0.001 

0.9 ±0.2 


All AGN, 2 4 3.66 

288 

1.27 ±0.05 

2 .8% 

0.552 ±0.017 

0.003 

1.1 ±0.3 


All AGN, 2 4 3.66, no LSS 

264 

1.31 ±0.05 

3.0% 

0.556 ±0.018 

0.003 

1.2 ±0.3 


RQ AGN 

172 

1.93 ± 0.11 

12 .8% 

0.709 ±0.022 

< 0.001 

2.3 ±0.1 


RQ AGN, 2 4 3.66 

155 

1.54 ±0.07 

5.2% 

0.708 ±0.023 

< 0.001 

3.0 ±0.2 


RQ AGN, 2 : ^ 3.66, no LSS 

141 

1.61 ± 0.08 

5.7% 

0.711 ±0.024 

< 0.001 

3.0 ±0.2 


RL AGN 

136 

1.00 ±0.06 

1.5% 

0.354 ± 0.025 

< 0.001 

-3.9 ±0.9 

-2.1 ± 0.3 

RL AGN, no LSS 

126 

1.03 ±0.07 

1 .6% 

0.362 ± 0.026 

< 0.001 

-3.6 ±0.9 

-2.0 ±0.3 


“Pure luminosity evolution P(z) = P(0)(1 ± z) kL 
6 Pure density evolution = <J>(0)(1 + z) kD 


14/14 distribution is similar to a uniform one in column 
(6), and the best fit parameters fez, and fco (when appli¬ 
cable) in column (7) and (8) respectively. When the best 
fit indicates negative luminosity evolution (i.e., fez, < 0), 
in fact, we fit a PDE model as well, which we feel is 
more physical in this case. The mean redshift is calcu¬ 
lated taking into account the effect of the sky coverage, 
that is each object is weighted by the inverse of the area 
accessible at the flux density of the source (e.g. Padovani 
et al. 2007). 

The fractional redshift distributions for the different 
classes are shown in Fig. 2. As also shown in Table 2, RQ 
AGN have the highest ( z) (z ~ 2) and reach z ~ 6 — 7 
(photometric redshifts). RL AGN have also a relatively 
broad distribution, extending up to z ~ 4.5 but with a 
lower (z) than RQ ones (z ~ 1). SFG are in between 
but closer to RL AGN. The main results on the sample 
evolution are the following: 

(i) The whole sample has (14/14) > 0.5 and shows 
significant evolution characterized by fez, = 1.6 ± 0.1; 

(ii) SFG evolve at a very high significance level (p- 
value < 0.001); their evolutionary parameter for the case 
of pure luminosity evolution is fen = 2.1 ± 0.1. This in¬ 
creases to fen = 2.4 ± 0.2 for a redshift cut z 4 3.25, 
which reduces the fraction of estimated redshifts. The 
evolution and LF of SFG will be presented and discussed 
in Padovani et al., in preparation; 

(iii) AGN as a whole evolve relatively weakly, with 
fen = 0.9 ± 0.2 (fen = 1.1 ± 0.3 for 2 4 3.66); 

(iv) RQ AGN, however, evolve very significantly (p- 
value < 0.001) with fen = 2.3 ± 0.1, consistent with the 
value of SFG (a thorough comparison between the two 
classes is deferred to Padovani et al., in preparation). This 
increases to fen = 3.0 ± 0.2 for a redshift cut z 4 3.66, 
which reduces the fraction of estimated redshifts; 

(v) RL AGN also evolve strongly (p-value < 0.001) 
but in the negative sense, with fen = —3.9 ± 0.9 and 
fcn> = —2.1±0.3. This evolution is in the opposite sense as 
that of RQ AGN. As we show below (Sect. 4.2 and 5.2.2), 
the best-fit evolution of RL AGN is more complicated 
than assumed here. 

We note that Padovani et al. (2011) found a dif¬ 
ference in (Ve/Va) between RL AGN with powers below 


and above P = 10 24 ’ 5 W Hz -1 , with the former hav¬ 
ing (Ve/Va) ~ 0.5 and the latter reaching smaller values. 
Given the flux density-limited nature of our sample, and 
the resulting degeneracy between power and redshift, we 
believe that this difference, which is still present in the 
current sample, is better characterised as a redshift de¬ 
pendence (see Sect. 5.2.2). 

Given the small area of our survey one could worry 
that the presence of large-scale structures (LSS) (e.g. Gilli 
et al. 2003) and related redshift spikes might influence 
some of our results. The LSS in the E-CDFS is discussed 
in Appendix B. To assess the impact of the LSS on our 
study, most of our results are presented for the “full” and 
the “no LSS” sub-samples, as illustrated, for example, in 
Fig. 2 (cf. the solid with the dotted lines). As shown in 
Table 2, the resulting (14/14) values and best-fit evolu¬ 
tionary parameters for the latter are well within lcr from 
those derived for the former, which shows that the ef¬ 
fect of the LSS on these results is minimal. These over¬ 
densities are obviously more noticeable when one studies 
the evolution of the LF with redshift (see, e.g. the top- 
middle and bottom-left panels of Fig. 9) but even then 
the revised LFs are typically within < 1 <t from the “full” 
ones. 


4.2 The banded (14/14) 

Given the available statistics, we can easily study pos¬ 
sible changes of the evolution with redshift by using 
the so-called banded (17/Vj„ax) statistic, i.e. the quan¬ 
tity ((V - 14)/(V4ax - 14)), or ((14 - Vo)/(Va - 14)) in 
our case, where 14 is the cosmological volume enclosed 
by a given redshift z 0 (Dunlop & Peacock 1990). This 
procedure allows us to detect any high-redshift (possibly 
negative) evolution by separating it by the well-known 
strong (positive) low-redshift evolution. 

Fig. 3 shows that ((14 - 14)/(14 - 14)) for RQ AGN 
has a very slow decline with redshift, but still with values 
significantly > 0.5, until it reaches 0.5 at z ft; 3.3. This 
suggests that the evolution might be somewhat weaker at 
higher redshifts; we will address this in Sects. 5.3.2 and 
6 . 1 . 2 . 

As for the RL AGN (14/14) is < 0.5 for the whole 
sample so the banded (V /Vjnax) statistic is not useful. 
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Figure 3. The banded statistic, ((14 — 14)/(14 — 14)) versus 
z 0 , for our RQ AGN with z 4 3.66. The horizontal dashed line 
indicates the value of 0.5, expected under the null hypothesis 
of no evolution. Open symbols take into account the presence 
of LSS as discussed in Appendix B. 

However, we can plot (14/14) vs. the maximum redshift 
of the sub-sample to see if there is evidence of positive 
or at least null evolution at low redshifts. This is done in 
Fig. 4, which shows that, indeed, up to z ~ 1.5 (14/14) 
is not significantly < 0.5, being actually > 0.5 at the 
~ 2 — 3cr level at z ~ 0.7 — 0.8. While part of this is due 
to the presence of LSS, taking these into account still 
results in (14/14) ~ 0.5 and actually > 0.5 at the 2. 60 - 
level at z ~ 0.7. At z > 1 (14/14) declines with redshift, 
reaching gradually its sample value. This shows that the 
RL AGN evolution changes from being positive/null to 
being negative around z ~ 1 — 1.5 and that therefore it 
cannot be characterized by a simple (l + z) fe law. In other 
words, (14/14) evaluated for the whole sample tells only 
part of the story. We will expand on this in Sect. 5.2.2. 

4.3 Maximum Likelihood Analysis 

A more general approach to estimate the evolution, and 
at the same time the LF, is to perform a maximum like¬ 
lihood fit of an evolving luminosity function to the ob¬ 
served distribution in luminosity and redshift. This ap¬ 
proach makes maximal use of the data and is free from 
arbitrary binning; however, unlike the 14/14 test, it is 
model dependent. We consider one and two power-law 
LFs, that is $(P) oc P' 71 and $(P) oc 1/[(P/P*) 71 + 
(P/P*) 72 ] respectively, with a break at P*. We refer the 
reader to Padovani et al. (2011) for more details. 

Our results are shown in Table 3, which gives the 


Figure 4. (14/14) as a function of maximum redshift for our 
RL AGN. The horizontal dashed line indicates the value of 
0.5, expected under the null hypothesis of no evolution, while 
the dotted line is the value for the sample. Open symbols take 
into account the presence of LSS as discussed in Appendix B. 

sample in column (1), the evolutionary model in column 
(2), the two LF slopes (if applicable) in columns (3) and 
(4), the best-fit evolutionary parameter in column (5), 
and the break power (if applicable) in column (6). Er¬ 
rors are la for one interesting parameter. The best-fit 
evolutionary parameters agree (within 1 — 2er) with those 
derived through the 14/14 approach. 

The maximum likelihood approach (but also the 
14/14 one) does not provide a goodness-of-fit test. There¬ 
fore, as it is typically done, we carry out a KS test on both 
the cumulative redshift and radio power distributions and 
compute the probability of having the KS test statistics 
as large or larger than the observed one. This can be used 
to reject a model when too low. The p-values are given 
in columns (7) (KS Z ) and (8) (KSp) respectively, where 
the > 2a and > 3<r levels correspond to values < 0.046 
and < 0.003 respectively. The fact that in most cases one 
of the two p-values is 4 0.04 shows that the models we 
adopted are too simplistic, as detailed in the next Section. 


5 LUMINOSITY FUNCTIONS 

The results of the maximum likelihood approach expand 
upon and complement those of the 14/14 test. A simple 
power law fit for the LF of AGN is excluded with very 
high significance (p -value < 0.001), both in the case of 
PLE and PDE. However, even a double power law fit is 
inconsistent with the data, which means that the adopted 
evolutionary laws are too simple. A simple power law 
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Table 3. Sample Luminosity Functions and Evolution: maximum likelihood analysis best fits. We consider a LF of the type 
4>(P) oc 1 /[(P/P *) 71 + {P/P*) 12 } and an evolution oc (1 ± z) k . KS Z and KSp are the KS p-values for the redshift and radio power 
distributions respectively. 


Sample 

Model 

71 

72 

k 

logP* 

W Hz ” 1 

KS Z 

KS P 

AGN, z < 3.66 

PDE 

1.70 ±0.04 


0.11 ± 0.26 


< 0.001 

0.016 

AGN, z ^ 3.66 

PLE 

1.70 ±0.04 


0 16+ 0 ’ 33 
U - 1O -0.41 


< 0.001 

0.016 

AGN, z ^ 3.66 

PDE 

±- — 0.2 

1.82 ±0.06 

0.0 ±0.3 

22.6 ± 0.6 

0.006 

0.306 

AGN, z < 3.66 

PLE 

1.2 ± 0.2 

1.80 ±0.05 

-0.15 ± 0.35 

22.6 ± 0.6 

0.005 

0.340 

AGN, 2 : ^ 3.66, no LSS 

PLE 

1.3 ± 0.2 

bo 

0 

1 + 

0 © 

b m 

cn 0 

0.1 ±0.4 

ZZ.O_Q 9 

0.039 

0.383 

RQ AGN, z ^ 3.66 

PLE 

2.3 ± 0.1 


2.8 ± 0.2 


0.022 

0.038 

RQ AGN, z 3.66 

PLE 

0 6 +0 ' 6 
u '°-0.7 

2 64”*”°' 19 
z -°^-0.27 

O C 9 + 0.20 

Z.OZ—0 23 

99 09 -I-O. 2 O 
ZZ.oz_q 1 e 

0.372 

0.885 

RQ AGN, z ^ 3.66, no LSS 

PLE 

0 9+ 0,5 
u ' —0.6 

9 7 Q + 0.25 

' z -0.23 

2.54 ±0.23 

22.50 ±0.25 

0.471 

0.879 

RL AGN 

PLE 

1 40 - *" 0-05 
J "* u — 0.04 

-6.0 ± 1.4 


0.003 

0.771 

RL AGN 

PDE 

1.40 ±0.04 


-2.4 ±0.3 


0.003 

0.772 

RL AGN, no LSS 

PDE 

1.41 ±0.04 


-2.3 ±0.3 


0.018 

0.870 


Models consistent with the data are indicated in boldface. 


fit for the LF of RQ AGN is also excluded (p-value ~ 
0.02), while a PLE evolution plus a two power-law LF 
is consistent with the data. As for RL AGN, neither a 
simple PLE or PDE model can fit the data, while a two 
power-law LF is not needed (the best fit slopes turn out 
to be the same well within the errors). The two AGN 
classes, apart from displaying evolutions with different 
signs, have also widely different LF shapes, with RL AGN 
being characterised by a much flatter LF. Note that the 
effect of the LSS on these results is truly minimal (< lcr). 

We now turn to study the AGN LFs in more detail. 
We first study all AGN together but then, given their 
widely different LFs and evolutions, we examine the two 
classes separately. 

The effect of cosmic variance was included in the LF 
uncertainties as follows: 1. we first estimated <Tdm for the 
relevant redshift range following Moster et al. (2011) (us¬ 
ing the IDL code provided in their paper); 2. we then 
derived the bias b for the appropriate redshift and num¬ 
ber density from Fig. 3 (left) of Somerville et al. (2004). 
We could not use in fact Moster et al. (2011) to get b as 
well given that their numbers apply to the whole galaxy 
population above a given stellar mass, while we here deal 
with a radio-selected sample; 3. we estimated the rela¬ 
tive cosmic variance as a v = b x <7dm (see Moster et al. 
2011, for details and definitions); 4. a v was then added in 
quadrature with the Poisson uncertainties for each sub¬ 
population in each redshift bin. We note that, since the 
number of objects per power bin in our LFs is never very 
large, the former turned out always to be smaller than 
the latter. 


5.1 All AGN 

Fig. 5 shows three different estimates of the local LF for 
our AGN, done using the 1/Vjnax, or in our case 1/14, 
method. The AGN z ^ 0.4 sample, which includes 29 
objects, is consistent with no evolution ((14/14) = 0.57± 
0.05) and can therefore be used as a robust proxy for the 
local AGN LF (filled squares; similar results are obtained 
by using the z ^ 0.3 sample, which however contains only 
21 objects). To increase our statistics, we have included in 


the sample all sources above our radio completeness limit 
independently of their f 3 . 611 m values (see Sect. 2.1), which 
adds four more sources 4 . In Fig. 5 we also plot the best 
fits to the local (z ^ 0.3) LFs from Sadler et al. (2002) 
and Mauch & Sadler (2007) (dotted and dash-dotted lines 
respectively). Our local LF is consistent within the errors 
(xt /S 1-5) with both but nominally a factor ~ 1.7 — 2.3 
higher than the Sadler et al. (2002) and Mauch & Sadler 
(2007) LF respectively. 

Our z < 0.4 LF agrees with the z < 0.2 AGN LF 
of Mao et al. (2012), based on the Australia Telescope 
Large Area Survey (ATLAS) at 1.4 GHz, which currently 
reaches ~ 150/rJy per beam. The Smolcic et al. (2009) 
AGN LF in the 0.1 — 0.35 redshift range (not shown for 
clarity) is also consistent with ours. 

Fig. 5 shows two further estimates of the local LF for 
E-CDFS AGN, which are based on the whole sample but 
are model dependent. The first one is the LF de-evolved 
to z = 0 (filled circles) using the best fit evolutionary pa¬ 
rameter from the (14/14) analysis. This LF is not consis¬ 
tent (xt ~ 6.1 — 8.2 respectively, p-value < 0.0001) with 
both previous LFs, which indicates that the redshift evo¬ 
lution is more complex than assumed. The situation does 
not change taking into account the LSS {xt ~ 5.0 — 6.2, 
p -value < 0.0001). The second one is the LF derived from 
the maximum likelihood analysis (double power-law fit, 
PLE), including also the LSS effect, which is ~ 2 — 3 
times above previous determinations. 

This situation is somewhat at variance with that in 
Padovani et al. (2011), where the local AGN LF was dif¬ 
ferent (p -value ~ 0.011) from both the Sadler et al. (2002) 
and Mauch & Sadler (2007) LF. We attribute this differ¬ 
ence to a statistical fluctuation due to the smaller (15 
sources vs. 29) local AGN sample available in that pa¬ 
per. As regards the fact that our LFs are still somewhat 
above previous LFs, we noted in Padovani et al. (2011) 
that our AGN LF includes the contribution of RQ AGN, 

4 We do not need in fact to restrict ourselves to sources with 
/ 3 . 6 /jm > 1 microjy to increase redshift completeness, since 
none of the AGN with z < 0.4 have redshift estimated from 

/3.6/im ■ 
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Figure 5. The local differential 1.4 GHz LF for AGN in a 
P X cj>{P) form. Filled (magenta) squares indicate the E-CDFS 
LF at z ^ 0.4, while filled (black) circles denote the LF de- 
evolved to z = 0 using the best fit evolutionary parameter from 
the (14/14) analysis, with open circles taking into account the 
presence of LSS. The solid line is the best-fit double power- 
law LF from the maximum likelihood method (PLE), while 
the dashed line is the same LF but taking into account the 
presence of LSS. Open (blue) squares refer to the z < 0.2 
ATLAS LF (Mao et al. 2012). The best fits to the local AGN 
LF from Sadler et al. (2002) (converted to our value of Hq) and 
Mauch & Sadler (2007) are also shown (blue dotted and dash- 
dotted lines respectively). Error bars for the E-CDFS z 4 0.4 
sample correspond to lcr Poisson errors evaluated using the 
number of sources per bin with redshift determination added 
in quadrature with the cosmic variance uncertainties. 

which were not present in significant numbers in these 
two LFs as these were based on the NRAO VLA Sky 
Survey (S'min > 2.8 mjy), while RQ AGN make up a 
non-negligible fraction of radio sources only below « 1 
mjy (Fig. 1). Moreover, both previous LFs included only 
non-stellar optical sources (i.e., they excluded quasars), 
while we have no such cut. The issue of the RQ AGN con¬ 
tribution can be dealt with by studying RL AGN only. 

5.2 RL AGN 

5.2.1 The local LF 

Fig. 6 shows three different estimates of the local LF 
for our RL AGN. The RL AGN z 4 0.4 sample, 
which includes 23 objects, is consistent with no evolu¬ 
tion ((14/14) = 0.57±0.06) and can therefore be used as 
a robust proxy for the local RL AGN LF (filled squares). 
As done for the full AGN sample, this includes all sources 
above our radio completeness limit independently of their 
/ 3 . 6 nm values (see Sect. 2.1). Our local LF is consistent 
within the errors (x% 1-3) with both the LFs from 

Sadler et al. (2002) and Mauch & Sadler (2007) (dotted 



Figure 6 . The local differential 1.4 GHz LF for RL AGN in 
a P X <f>(P) form. Filled (magenta) squares indicate the E- 
CDFS LF at z 4 0.4, while filled (black) circles denote the LF 
de-evolved to z = 0 using the best fit evolutionary parameter 
from the (14/14) analysis, with open circles taking into ac¬ 
count the presence of LSS. The solid line is the best-fit from 
the maximum likelihood method (ZDE), while the dashed line 
is the same LF but taking into account the presence of LSS. 
Open (blue) squares refer to the 0.05 < z < 0.35 VIDEO- 
XXM3 LF (McAlpine, Jarvis, & Bonfield 2013), while open 
(red) triangles are the z 4 0.5 Subaru /XMM-Newton Deep 
Field LF (Simpson et al. 2012), corrected to take into account 
the different selection criteria as described in the text. The 
best fits to the local AGN LF from Sadler et al. (2002) (con¬ 
verted to our value of Hq), Mauch & Sadler (2007) (blue dotted 
and dash-dotted lines respectively) and Best et al. (2014) (all 
sources: red short-long dashed line; radiative-mode: green dot 
- short-dashed line) are also shown. Error bars for the E-CDFS 
z 4 0.4 sample correspond to lcr Poisson errors evaluated us¬ 
ing the number of sources per bin with redshift determination 
added in quadrature with the cosmic variance uncertainties. 


and dash-dotted lines respectively), although still a factor 
~ 1.4 — 1.8 higher. 

We also compare our local LF with the z 4 0.5 
Subaru /XMM-Newton Deep Field LF of Simpson et al. 
(2012) (red open triangles). That paper adopts a some¬ 
what different method to classify AGN into RL and RQ 
sources, based on the rest-frame ratio of mid-IR to ra¬ 
dio flux density (724 = log(S24 M m/5i.4GHz), with RL AGN 
defined by 524 4 —0.23. Our method relies instead on 
the observed 924 and the dividing line is redshift depen¬ 
dent. To take this difference into account, we checked how 
many RL AGN with z 4 0.5 and flux density ^ 100 ^tJy 
(the limit of Simpson et al. 2012) were classified as RL 
AGN by our method and by that of Simpson et al. (2012) 
and corrected their LF by the ratio of the two (1.6). 
After doing this, the two LFs are in perfect agreement. 
Note that Del Moro et al. (2013) have also shown that 
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Table 4. RL AGN Luminosity Function and Evolution: ZDE maximum likelihood analysis best fits. We consider a LF of the type 
4>(P) oc P —71 and an evolution peaking at z pea k, which is oc (1 + z) k i°w and oc (1 + z) khi s h at low and high redshifts respectively 
(see eq. 1). KS Z and KSp are the KS p- values for the redshift and radio power distributions respectively. 

Sample 'yi ^-low ^-high -^peak KSz KSp 

RL AGN 1.44 ±0.04 2.2±^;| -3.9±°'g 0.49 ±0.10 0.331 0.889 

RL AGN, no LSS 1.45 ±0.05 1.4 ° -3.715;! 0.41 ±0.18 0.524 0.872 


524 ^ —0.23 gives an incomplete selection of RL AGN as 
it misses evidence for radio jets. 

We also include the 0.05 < z < 0.35 VIDEO -XXM3 
LF of McAlpine, Jarvis, & Bonfield (2013) (blue open 
squares), who selected their AGN based on them being 
in host galaxies “redder than the spiral galaxy templates” 
(and therefore are most likely going to be of the RL type). 
Their LF appears to be somewhat above ours at Pm 10 23 
W Hz -1 but the difference in selection criteria between 
the two derivations makes a more detailed comparison 
difficult. 

Finally, Fig. 6 shows also two more LFs: the z < 0.3 
LF of Best et al. (2014) for all sources (red short-long 
dashed line) and radiative-mode only (green dot - short- 
dashed line). Our local LF is consistent with the first 
LF although nominally a factor ~ 1.3 times higher. We 
discuss this, together with the general issue of comparing 
our LFs to many previous derivations, in Sect. 6.1.1. Our 
LF is ~ 2 orders of magnitude above the radiative-mode 
radio AGN LF, which confirms the fact that the bulk of 
our RL AGN are of the jet-mode type (see also Sect. 6.2). 

Fig. 6 shows also two further estimates of the local 
LF for E-CDFS RL AGN, which are based on the whole 
sample but are model dependent. The first one is the 
LF de-evolved to z = 0 (filled circles) using the best fit 
evolutionary parameter for PDE from the (14/14) anal¬ 
ysis. This LF is totally inconsistent (p-value < 0.0001) 
with both previous LFs, which indicates that the red- 
shift evolution is more complex than assumed. (The situa¬ 
tion does not change taking into account the LSS: p -value 
< 0.0001). This is confirmed by the second one, that is 
the LF derived from the maximum likelihood analysis us¬ 
ing a more complex density evolutionary law than PDE, 
as detailed in Sect. 5.2.2. In this case the best fit LF is in 
good agreement with both previous LFs up to P > 10 24 
W Hz -1 , which suggests that this evolutionary law might 
be valid only up to these powers, as we cannot yet con¬ 
strain the bright end slope of the LF (this has however 
been done by other wider/shallower surveys). 

5.2.2 The LF redshift evolution 

We noticed in Sect. 4.2 that the RL AGN evolution can¬ 
not be simply parametrized by a (l±z) / ' law but reaches 
a peak at z ~ 0.5— 1. We therefore introduce a more com¬ 
plex density evolutionary law, which we term as “Zpeak 
density evolution” (ZDE) as follows: 

(1 + Zc)* 1 ™ + (1 + Zc)] fchigh ,, 

[(1 ± 2 C )/(1 ± «)] fc low ± [(1 ± z c )/(l + 2 )] fc high ' 

This allows for a situation where evolution changes 
sign above some redshift and for a smooth transition. In- 



22 24 26 22 24 26 22 24 26 


log P r (1.4 GHz) [W Hz -1 ] 

Figure 7 . The differential 1.4 GHz LF for RL AGN in a 
P x (f)(P ) form in six redshift bins. The (blue) solid lines repre¬ 
sent the best fit single power-law LF from the maximum likeli¬ 
hood method (ZDE) evolved to the central redshift of the bin 
using the best fit for the complex density evolution described 
above, with dotted lines showing the same LF at the two ex¬ 
treme redshifts defining the bin. The (red) short-dashed line 
represents the best fit LF at 2 = 0. Error bars correspond to 
lcr Poisson errors evaluated using the number of sources per 
bin with redshift determination added in quadrature with the 
cosmic variance uncertainties. The percentage of redshifts es¬ 
timated from /3.6/i,m is also given for each bin. Open symbols 
in the z = 0.5 — 1.0 bin take into account the presence of LSS, 
while (green) triangles are the LF from Best et al. (2014) (all 
sources). The dotted vertical lines denote the completeness 
limit in power for each redshift range. See text for details. 


deed, ev(z) oc (1 + z) fclow and ev(z) oc (1 + z ) khish for 
z <C z c and z z c respectively (and eu(0) = 1). If 
faow > 0 and khigh < 0 then ev(z) peaks at z pe ak = 
(—fehighAiow) 1 ^^ fchigh ’” fclow ^(l + z c ) — 1. The results of fit¬ 
ting this model to our RL AGN sample are reported in 
Table 4. We find a best fit z pea k ~ 0.5, with the density 
evolution turning from positive to negative respectively 
below and above this value, while the slope of the LF 
is basically the same as derived before. Importantly, the 
model is consistent with the data as shown by the results 
of both KS tests. Taking into account the effects of the 
LSS changes our best fit parameters by less than la. 
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Table 5. Luminosity Functions of VLA-ECDFS RL AGN. 
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Numbers in parenthesis take into account the effect of the 
LSS. Errors correspond to 1 rj Poisson errors evaluated using 
the number of sources per bin with redshift determination 
added in quadrature with the cosmic variance uncertainties. 
The conversion to units of Mpc -3 dex- 1 used, for example, by 
Simpson et al. (2012), is done by subtracting 9 — log(ln(10)) 
from our values. 


Fig. 7 (tabulated in Table 5) shows the evolution of 
the RL AGN LF in six redshift bins. The decrease in 
number at higher redshifts is clearly visible, as is the fact 
that the ZDE maximum likelihood fit provides a very 
good representation of the LF redshift evolution. 

As discussed in Sect. 6.1.1, a direct comparison be¬ 
tween the evolution of our LF and previously published 
results is not straightforward. However, once one takes 
into account some of the most obvious differences our re¬ 
sults are in agreement with previous ones. As an example, 
Fig. 7 shows that our LF in the 0.5 — 1.0 redshift bin is 
somewhat higher, as expected, but still consistent within 
the errors, with that of Best et al. (2014) (all sources). 


The same thing applies to the LFs of Gendre, Best & 
Wall (2010) in the 0.8 — 1.5 and 1.2 — 2.5 bins. 

5.3 RQ AGN 

5.3.1 The local LF 

The local (z < 0.3 — 0.4) LF for radio-selected RQ AGN 
has, so far, never been determined 5 . We present it in Fig. 
8. The RQ AGN z ^ 0.4 sample (filled squares) includes 
only 6 objects and the corresponding LF is therefore quite 
uncertain. We therefore complement it by deriving the 
LF for the full sample de-evolved to 2 = 0 (filled cir¬ 
cles) using the best fit evolutionary parameter from the 
(Ve/Va) analysis and the local LF derived from the max¬ 
imum likelihood analysis (2 power-law, PLE). Given the 
discussion in Sect. 4 both these estimates should be reli¬ 
able. We also plot in Fig. 8 the radio LFs for the 12/rm 
(open circles) and Center for Astrophysics (CfA; open 
squares) Seyferts (Rush, Malkan, & Edelson 1996) and 
the Palomar Seyferts (Ulvestad & Ho 2001), that is of 
the local counterparts of (mostly) RQ AGN. There are a 
couple of points one needs to consider here: 1) all these 
LFs are bivariate, that is they were derived by observing 
in the radio band samples of Seyfert galaxies. Given that 
the original sample selection was done in the optical, it 
is almost certain that the samples are not complete in 
the radio and therefore the LFs are lower limits to the 
true ones; 2) there is some scatter between the various es¬ 
timates, although this is usually within the rather large 
error bars, which are due to the fact that many bins in¬ 
clude only 1 — 2 sources; 3) the high-power (P > 10 24 
W Hz -1 ) tail should not be considered since the 12/rm 
sample, for example, includes a number of “classical” RL 
sources (e.g. 3C 273, 3C 120, OJ 287). Taking all of the 
above into account, our determination of the local LF 
for radio-selected RQ AGN is not inconsistent with the 
Seyfert LFs. We will return to this issue in Sect. 6.1.2. 

5.3.2 The LF redshift evolution 

The first ever determined evolution of the RQ AGN radio 
LF is shown in six redshift bins in Fig. 9 (tabulated in Ta¬ 
ble 6). The increase in power at higher redshifts is clearly 
visible, as is the fact that the PLE maximum likelihood 
fit provides a good representation of the LF redshift evo¬ 
lution, which means that the possible slowdown at higher 
redshifts hinted at in Sect. 4.2 has to be subtle. 

5.4 RL and RQ AGN comparison 

Fig. 10 shows the redshift evolution of RL and RQ AGN 
in various luminosity bins, with the points coming from 
the 1/14 method and the lines being the best maxi¬ 
mum likelihood fits assuming PLE for RQ AGN and 
PDE for RL AGN respectively and a single power law 
(given the relatively small power range). The figure shows 
that, while at low powers (Pi ,4GHz < 10 23 W Hz" 1 ) 

5 In Padovani et al. (2011) we estimated it by de-evolving the 
full LF to z = 0. 
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Figure 8. The local differential 1.4 GHz LF for RQ AGN in 
a P X (pi P) form. Filled (magenta) squares triangles indicate 
the E-CDFS LF at z ^ 0.4, while filled (black) circles denote 
the LF de-evolved to z = 0 using the best fit evolutionary 
parameter from the {Ve/Vs) analysis. The solid line is the best- 
fit double power-law LF from the maximum likelihood method 
(PLE), while the dashed line is the same LF but taking into 
account the presence of LSS. The (green) dotted line is the 
best fit to the local LF of radiative-mode radio AGN (Best et 
al. 2014). Local bivariate LFs for Seyfert galaxies from three 
samples are also shown: the 12pm (open blue circles) and CfA 
(open red squares) Seyferts (Rush, Malkan, & Edelson 1996) 
and the Palomar Seyferts (filled cyan triangles; Ulvestad & Ho 
2001). Error bars for the E-CDFS z ^ 0.4 sample correspond 
to ler Poisson errors evaluated using the number of sources 
per bin with redshift determination added in quadrature with 
the cosmic variance uncertainties. 

RQ and RL AGN are equally numerous, the situation 
changes at 10 23 < Pi.4GHz < 10 25 W Hz” 1 , where RQ 
AGN overtake RL ones. For example, at z ~ 2 — 3 and 
10 24 < Pi.4GHz < 10 25 W Hz -1 RQ AGN are more abun¬ 
dant than RL ones by factors ~ 3 — 10. At even larger 
powers (Pi.4GHz > 10 25 W Hz -1 ) the number densities 
of the two classes appear to be similar again but at this 
point we are running out of RQ AGN. 

We stress that the fact that the decline appears to 
start at larger redshifts for more powerful sources has 
nothing to do with “cosmic down-sizing” (as instead 
stated by Simpson et al. 2012) but is simply an effect 
of our flux density limit. Within a redshift bin, in fact, 
not all sources in a given radio luminosity bin might be 
visible by our survey and this leads to a “loss” of ob¬ 
jects. As a result, sources will start to fall below the flux 
density limit at progressively higher redshifts for larger 
powers. This explains the sudden stop in the rise of both 
classes, which in deeper samples is expected to continue 
(obviously barring an intrinsic redshift cutoff). 

Ueda et al. (2014) reach maximum number densities 
~ 3 x 10 5 Gpc” 3 for their X-ray selected (largely RQ) 


Table 6. Luminosity Functions of VLA-ECDFS RQ AGN. 


z range l°gPi.4GHz logP4>(P) N 

W Hz” 1 Gpc” 3 




21.77 

k no+0.37 
°- u — 0.49 


2 

z 0.5 


22.17 

4 24"*" 0 ' 52 


1 



22.57 

4 43+ 0 - 24 
0.31 


5 



22.97 

4 21"*" 0,37 
0.53 


3 



22.88 

4 85”*” 0 ' 12 

^• oo -0.14 

(4.48±°+ 9 ) 

18 (10) 

O 

Cn 

A 

N 

O 

V/ 

23.28 

4 46 -1 " 0,12 

-0.14 

(4.43+gJ 2 ) 

20 (19) 



24.08 

3 09”*” 0,52 
0.91 

1 



24.48 

3 09"*~ 0 ’ 52 

o.u»_ 0 91 


1 



23.50 

4 65 +0 ' 13 

^• uo -0.15 


18 

o 

A 

N 

s; 1.5 

23.90 

^.uo_o is 


14 



24.30 

3 59"*" 0 ’ 26 

o.D»_0 33 


5 



23.57 

q cq+0.37 
o.oo_ o 49 

(4-03±°+) 

2 

1.5 < z 

o 

CN 

V/ 

23.97 

4 06 +0 ' 14 
^• uo -0.16 

14 (13) 



24.37 

q '71 +0.18 
' -0.22 

(3.66±°' 2 °) 

8 (7) 



24.77 

q 11 +0.37 

( 2 .81±g:|2) 

2 (1) 



24.29 

3 97 - *" 0,16 

1 -0.19 


11 

2.0 < z 

lO 

CN 

V/ 

24.69 

o.ui7_o 52 


2 



25.09 

9 70+0.52 

0.91 


1 



25.49 

9 7Q+0.52 

1 0.91 


1 



24.32 

q 7C+O.I8 
<o -0.20 


9 

2.5 < z 

< 3.66 

24.72 

3 «Q + 0 - 15 

o.°o_0 17 


14 



25.12 

9 7q+0.37 
0.48 


2 



25.52 

2 43+ 0 52 
0.83 


1 


Numbers in parenthesis take into account the effect of the 
LSS. Errors correspond to la Poisson errors evaluated using 
the number of sources per bin with redshift determination 
added in quadrature with the cosmic variance uncertainties. 
The conversion to units of Mpc” 3 dex” 1 used, for example, by 
Simpson et al. (2012), is done by subtracting 9 — log(ln(10)) 
from our values. 


AGN with 10 42 < I/ x < 10 43 erg s~ 4 and z ~ 0.7 (see 
their Fig. 12), to be compared with our value of ~ 7 x 10 4 
Gpc -3 at z ~ 0.3 (Fig. 10). This is to be expected as their 
AGN surface density is ~ 3, 000 deg -2 vs. ~ 860 deg” 2 
for our RQ AGN (see Sect. 6.1.2), which simply means 
that we have not detected all AGN with L x > 10 42 erg 
s' 1 (but only those above our radio flux density limit). 

Finally, we note that the number densities for “clas¬ 
sical” RL quasars are orders of magnitude smaller than 
those of our RL AGN (< 1 Gpc" 3 for P 2 . 7 ghz > 10 27 W 
Hz” 1 : see, e.g. Fig. 5 of Wall et al. 2005). 


6 DISCUSSION 

6.1 The evolution of micro-Jy radio sources 

We now analyze in more detail the evolution of the two 
AGN classes and, when applicable, compare it with pre¬ 
vious results. 
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log P r (1.4 GHz) [W Hz" 1 ] 

Figure 9. The differential 1.4 GHz LF for RQ AGN in a 
P X 0( P) form in six redshift bins. The (blue) solid lines rep¬ 
resent the best fit double power-law LF from the maximum 
likelihood method evolved to the central redshift of the bin 
using the best fit for pure luminosity evolution (1 + z) 2 ' 52 , 
with dotted lines showing the same LF at the two extreme 
redshifts defining the bin. The (red) short-dashed line rep¬ 
resents the best fit LF at z = 0. Error bars correspond to 
lcr Poisson errors evaluated using the number of sources per 
bin with redshift determination added in quadrature with the 
cosmic variance uncertainties. The percentage of redshifts es¬ 
timated from /3.6(im is also given for each bin. Open symbols 
in the z = 0.5 — 1.0 and 1.5 — 2.0 bins take into account the 
presence of LSS. The (black) dotted vertical lines denote the 
completeness limit in power for each redshift range. See text 
for details. 


6.1.1 RL AGN 

The evolution of RL AGN is not described by PLE or a 
simple PDE. Instead, it can be characterised as a density 
evolution reaching a peak at z ~ 0.5 and then decreasing. 

Rigby et al. (2011) studied various samples of steep- 
spectrum (a r > 0.5) AGN selected from a variety of ra¬ 
dio surveys with increasingly smaller areas and flux den¬ 
sity limits. They found a significant decline in comov¬ 
ing density at z > 0.7 for their lower luminosity sources 
(10 25 < Pi.4GHz < 10 26 W Hz -1 ). Turnovers were still 
present at higher powers but moved to larger redshifts, 
suggesting a luminosity-dependent evolution of the red¬ 
shift peak, similar to the so-called “cosmic downsizing” 
seen, for example, in the X-ray band (Ueda et al. 2014). 
Fig. 10 (bottom right) shows the relevant points (filled 
triangles) from Fig. 11 of Rigby et al. (2011) in one of the 
two ranges of luminosities we have in common. Consider¬ 
ing that, unlike these authors, we have no cut on a r and 
therefore we are bound to have more AGN in our sample, 
our points agree with theirs within our error bars. The 
same is true for the 10 26 —10 27 W Hz -1 range (not shown 
because there are no RQ AGN at those powers). Fig. 10 



0 1 2 3 0 1 2 3 

Z 


Figure 10. The redshift evolution of RL (red squares) and 
RQ (blue circles) AGN in various luminosity bins. The top 
labels refer to cosmic time in Gyr. The points come from the 
1 /Va method, while the lines are from the best fits assuming 
a single power law and PLE for RQ AGN and PDE for RL 
AGN respectively, taking into account the effect of our flux 
density limit. Error bars correspond to lcr Poisson errors added 
in quadrature with the cosmic variance uncertainties. Empty 
symbols refer to the “no LSS” samples. Green triangles are 
from Rigby et al. (2011). 

shows that at the powers corresponding to the low lumi¬ 
nosity end of the Rigby et al. (2011) sample we cannot 
put any constraints on a peak at z ~ 0.7, as we are only 
sensitive to z > 0.5, but our data are consistent with the 
trend present in the Rigby et al. (2011) sample. Moreover, 
we can use Fig. 10 of Rigby et al. (2011), which shows 
a strong correlation between z pea k and power. Since our 
median Pi.4GHz is ~ 10 24 W Hz -1 the z pea k ~ 0.5 we find 
through the ZDE fit appears to be in agreement with the 
Rigby et al. (2011) results. 

Best et al. (2014) have also studied various samples 
of (steep-spectrum) radio AGN selected from a variety 
of surveys with increasingly smaller areas and flux den¬ 
sity limits. They classify their sources into radiative-mode 
and jet-mode AGN using emission line diagnostics. Based 
on Fig. 6, Bonzini et al. (2013), and Sect. 6.2 (see also 
Fig. 1 of Best et al. 2014), the large majority of our 
RL AGN have to be of the latter type, so we compare 
our findings to their results on the jet-mode population. 
For Pi.4GHz < 10 24 W Hz 1 Best et al. (2014) (see their 
Fig. 5) find that its space density stays constant up to 
z ~ 0.5 and then decreases to z = 1. At moderate pow¬ 
ers, 10 24 < Pi. 4 GHz < 10 26 W Hz -1 , the space density 
increases to z ~ 0.5 before falling. At the highest powers 
the space density appears to increase up to z ~ 1 but the 
statistics is somewhat limited. 

Our results are similar to those of Best et al. (2014), 
taking into account our somewhat more limited coverage 


© 2015 RAS, MNRAS 000, 1-19 




































14 P. Padovani et al. 


of the luminosity - redshift plane. Fig. 10 shows that for 
Fi. 4 ghz < 10 23 W Hz -1 the space density of RL AGN in¬ 
creases up to ^ K 0.5 while for 10 24 < Pi.4GHz < 10 26 W 
Hz -1 , and taking into account the effect of the LSS, it in¬ 
creases to z ~ 1.0 before falling. We cannot say anything 
about larger powers. 

It is important to stress the difference between a 
“purely flux density limited” sample, like ours, and many 
others published in the literature, which are very often 
put together by cross-correlating very large surveys. For 
example, the Best et al. (2014) LF, as discussed in Heck¬ 
man & Best (2014), is built by combining four different 
radio LFs based on four spectroscopic samples: Las Cam- 
panas (Machalski & Godlowski 2000), 2dFGRS (Sadler et 
al. 2002), SDSS/NVSS/FIRST (Best & Heckman 2012), 
and 6 dFGS (Mauch & Sadler 2007). Although these sam¬ 
ples are large (including, e.g. more than 2,500 and 7,000 
sources for the 6 dFGS and SDSS/NVSS/FIRST sam¬ 
ples) one needs to consider that, by construction, they do 
not sample the full RL AGN population. All of them, in 
fact, are limited to galaxies (and therefore do not include 
non-stellar and broad-line sources) and to steep-spectrum 
(a r > 0.5) radio sources. Gendre, Best & Wall (2010) and 
Gendre et al. (2013), who derived the LF for FRI/II radio 
galaxies, also included only galaxies, excluding any com¬ 
pact or quasar-dominated source. Our LF, being simply 
radio flux density limited, has no such biases and there¬ 
fore is bound to find number densities larger than these 
determinations; this is true also for Mao et al. (2012) and 
Simpson et al. (2012) (pending their different definition 
of RL AGN). 

Padovani et al. (2011) made a connection between 
the negative evolution of RL AGN and that of elliptical 
galaxies. For example, Taylor et al. (2009) have found 
that the number density of massive (M* > lO n M 0 ) 
red galaxies declines with redshift as 4?(.z) oc (1 + 

3)-1.60±0.14(±0.21) for ^ ^ 1 g6 j f wg rggtrict our gam _ 

pie to z < 1.8 we get <J>(^) oc (1 + ,z) -1 ' 3±0 ' 5 , which is 
very similar. 

Fig. 11 compares the LF of our RL AGN in the red¬ 
shift range 1.3 — 1.7 with two different models at 2 = 1.5 
from Best et al. (2014). These models, starting from the 
assumption that jet-mode RL AGN are hosted in quies¬ 
cent galaxies, combine the known stellar mass function 
of the host galaxies with the stellar mass dependency of 
the presence of a RL AGN to make predictions on the 
evolution of the RL AGN LF. Model 2b is a luminosity- 
density evolution model, in which the luminosity of the 
RL AGN increases with increasing redshift and there is 
a delay of ~ 2 Gyr between the formation of the host 
galaxy and the presence of a radio jet. Model 3a is a 
pure density evolution model in which a subset of the 
sources classified as jet-mode are in reality related to the 
radiative-mode RL AGN population. This could be due 
to a misclassification of some objects or to a change in 
physical properties with redshift. As shown in Fig. 7 of 
Best et al. (2014), the LFs predicted by these models are 
clearly distinguishable only at 2 > 1 , i.e. beyond their 


6 Best et al. (2014) have compiled some more recent results, 
which paint a somewhat more complex evolution. 
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Figure 11. The RL AGN LF at 1.3 C ~ 4 1.7 compared 
with two different models from Best et al. (2014) at z = 1.5. 
Filled (black) circles include all sources, while (red) triangles 
represent our best estimate of the jet-mode LF. See text for 
details. 

maximum redshift. Our sample covers a wide radio power 
range also at z ~ 1.5 and therefore we have the oppor¬ 
tunity to test these different models. As shown in Fig. 
11, model 3a is consistent with our LF, while model 2b 
is clearly not. Since the predictions were based on jet¬ 
mode only RL AGN, we also restricted our sample only 
to this type of sources by following Bonzini et al. (2013) 
and using the 22 /im power criterion proposed by Giirkan, 
Hardcastle, & Jarvis (2014) 7 . Only one of our jet-mode 
RL AGN turned out to have an Eddington ratio incon¬ 
sistent with this classification (> 0.1; see Sect. 6.2) and 
was therefore re-classified. The jet-mode LF is shown in 
Fig. 11 (filled red triangles) and is basically the same as 
the overall LF. 

6.1.2 RQ AGN 

We have provided a more robust estimate of the evolution 
of RQ AGN in the radio band, following upon the first one 
derived in Padovani et al. (2011). We have modelled it as 
a pure luminosity evolution P(z) oc (1 + z ) kL , obtaining 
fer, = 3.0±0.2 from the 14/14 analysis and fcz, = 2.5±0.2 
from the maximum likelihood one, in both cases over the 
0.2 — 3.66 redshift range. 

The redshift evolution of X-ray selected (largely RQ) 
AGN is best described by a more complex evolutionary 
model, the so-called luminosity-dependent density evo¬ 
lution (LDDE), which takes into account the observed 

7 As discussed in Bonzini et al. (2013), this method has not 
been tested on as faint a sample as ours and therefore we 
regard this as only a rough classification. 
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steepening with redshift in the faint-end slope of the LF 
and the fact that luminosity evolution appears to stop at 
z « 2.5, as shown, for example, in Fig. 11 of Ueda et al. 
(2014). It is important to note that the latter LF, and all 
those studied in the X-ray band, refer only to Compton 
thin AGN, characterised by an absorbing column density 
Nu < 1.5 x 10 24 cm -2 , while in principle we have no such 
limitation (although our cut in X-ray power might also 
exclude Compton thick sources). Keeping that in mind, 
Fig. 9 shows that the LF faint-end slope is constrained 
only at z < 0.5; moreover, most (~ 74%) of our RQ AGN 
have z < 2. The two effects combined make any devia¬ 
tion from PLE of the type seen in the X-ray band hard 
to detect. 

We nevertheless split the RQ AGN sample in two 
sub-samples of the same size with redshift respectively 
smaller and larger than 1.3, fitting PLE evolution to 
both. Evolution slows down at higher redshifts, with 
k = 4.0 ± 0.6 and k = 2.0 ± 0.5 respectively below and 
above z = 1.3. However, the errors on the evolutionary 
parameters overlap at the ~ 2cr level, making this differ¬ 
ence not statistically significant. In other words, a high- 
redshift slowing down of the evolution is consistent with, 
but not required by, our data. 

We note that the surface density of our RQ AGN, 
~ 860 deg -2 , is already larger than that reached by the 
deepest optical surveys, that is ~ 670 deg -2 for type I 
and II combined. This number is the sum of ~ 200 deg -2 , 
from a type II AGN sample based on the zCOSMOS sur¬ 
vey and selected on the basis of the optical emission line 
ratios (Bongiorno et al. 2010), and ~ 470 deg -2 , from 
a type I AGN sample extracted from the VIMOS VLT 
Deep Survey (Gavignaud et al. 2006). We note that the 
overall (RQ plus RL) AGN surface density is ~ 1, 500 
deg -2 . 

As regards the X-ray band, our RQ AGN surface 
density is ~ 20 times smaller than that of Lehmer et al. 
(2012), which is simply due to the fact that our radio flux 
density limit is still not as deep as the equivalent one in 
the X-rays (see also Fig. 10 in Vattakunnel et al. 2012, 
for a comparison of the star formation rates reachable in 
the two bands). 

To allow for a fairer comparison 8 between our sample 
and the X-ray samples we have fitted a double power 
law LF to all points in Fig. 8, obtaining Nt ~ 2.7 x 
10 5 Gpc -3 , a value, which is basically the same as the 
maximum number density of Ueda et al. (2014). This 
implies a surface density ~ 15, 400 deg -2 (using eq. 2 of 
Padovani 2011), to be compared with ~ 15, 000 deg -2 in 
the 4 Ms survey of Lehmer et al. (2012). In other words, 
once one takes into account the fact that our radio flux 
density limit is still not as deep as the equivalent one 
in the X-ray band by including the radio LF of Seyfert 
galaxies, the number and surface densities one obtains 
are the same, which shows that we are simply looking 
at the same population from two different bands and the 


8 Some of the topics discussed here and in Sect. 6.3 were orig¬ 
inally presented in Padovani (2011). We give here updated 
numbers. 
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Figure 12. The distribution of the ratio of X-ray power to 
stellar mass, L x /M+ for RQ (blue solid histogram) and RL (red 
dashed histogram) AGN. The top labels give the approximate 
values of the Eddington ratio. Shaded areas are upper limits. 

sources we are selecting in the radio band are the same 
as the X-ray emitting AGN. 

6.2 Astrophysics of micro-Jy sources 

To analyse in more depth the astrophysical nature of our 
AGN, we study the ratio of X-ray power to stellar mass, 
L x /M*, a proxy for the Eddington ratio L/L B dd (see, e.g. 
Bongiorno et al. 2012). This can in fact be written as 

L _ _ kboi _ L x . . 

L b dd ~ 1.26 x 10 38 T MR y ’ 

where fcbol is the 2—10 keV bolometric correction and 
T = Mbh/AZ* is the ratio between black hole and stellar 
mass. 

Stellar masses were estimated by Bonzini et al. 
(2013), using the method of Bongiorno et al. (2012). Fig. 
12 displays the distribution of L x /M* for RQ (blue solid 
histogram) and RL (red dashed histogram) AGN, with 
shaded areas denoting upper limits. We assume fcbol = 25 
(Bongiorno et al. 2012) and T = 0.001 (Trump et al. 
2015, and references therein), noting that, given the ob¬ 
servational scatter and a possible dependence of fcbol on 
L/L b ad, the resulting L/L B dd values are only approxi¬ 
mate. These are given by the top labels in Fig. 12. For 
our choice of parameters, L x /M* = 5 x 10 33 erg s -1 Mg 1 
corresponds to the Eddington limit. 

Fig. 12 shows that the two L x /M* distributions ap¬ 
pear very different, with RQ AGN having much larger 
values. Given the substantial fraction of upper limits (all 
on L x ) we used ASURV (Lavalley, Isobe, & Feigelson 
1992), the Survival Analysis package which employs the 
routines described in Feigelson & Nelson (1985) and Isobe 
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et al. (1986), which evaluate mean values by dealing prop¬ 
erly with non-detections and also compute the probabil¬ 
ity that two samples are drawn from the same parent 
population. RQ and RL AGN have totally different ( p- 
values < 0.001) distributions, according to all two-sample 
tests in ASURV. The means are also very different, being 
Lx/M* ~ 2.4 x 10 32 erg s _1 Mg 1 and ~ 5.5 x 10 3 ° erg 
s' 1 Mq 1 for RQ and RL AGN respectively, a factor ~ 40 
difference. The corresponding (approximate) Eddington 
ratios are ~ 0.05 and « 0.001. Given that the dividing 
line between radiative- and jet-mode sources is thought 
to be L/L Bdd ~ 0.01 (e.g. Heckman & Best 2014), our 
results are consistent with our RQ and RL AGN being 
mostly of the former and latter type, respectively. 

6.3 The intrinsic fraction of radiative-mode RL 

AGN 

The RL AGN fraction, /rlagn, has been derived in the 
past from bright, optically selected samples which contain 
RL sources very different from most of the sources studied 
here. Namely, they are RL quasars and therefore belong 
with the radiative-mode AGN, while we have mostly jet¬ 
mode objects (Sect. 6.2). Indeed, the RL AGN in the 
“classic” Palomar Green sample studied by Kellermann 
et al. (1989), from which the “standard” value /rlagn ~ 
10% comes, 9 have Psghz > 10 25 W Hz -1 and maximum 
number densities P<J?(P) ~ 11 Gpc -3 (Padovani 1993), 
that is they are below the plotting area in Figs. 5 and 6 
and barely inside in Fig. 8. 

The way this fraction was derived was by looking for 
a bimodality in the radio-to-optical flux density ratio, 
R, classifying as RL the sources with R larger than the 
value corresponding to the minimum in the distribution 
(0.1 < R < 100 in Kellermann et al. 1989) 10 . Padovani 
(2011) has pointed out that: a) bright, optically selected 
samples might include RL quasars with their optical flux 
boosted by relativistic beaming (e.g. Goldschmidt et al. 
1999), which artificially increases their fraction; b) RL 
AGN are on average more powerful than RQ ones in 
the optical band (e.g. Zamfir, Sulentic & Marziani 2008). 
When one reaches the very faint end of the optical LF 
only RQ AGN will be present and therefore the inte¬ 
grated, intrinsic RL fraction will be quite small. Both 
of these effects lead to an overestimate of the RL AGN 
fraction when using the “standard” method. Indeed, it 
has been known for some time that the radio-loud frac¬ 
tion drops with decreasing optical luminosity (Padovani 
1993). Jiang et al. (2007) have also shown that the frac¬ 
tion of RL quasars at z = 0.5 declines from 24% to 6% 
as luminosity decreases from M 2500 = —26 to —22. In 
other words, the usually quoted value refers to the bright 

9 Although in that paper the quoted value is 15 — 20%. 

19 We obviously cannot study the R distribution to look for 
a possible bi-modality as our RL AGN are mostly jet-mode 
sources while our RQ AGN are mostly radiative-mode sources 
(see Sect. 6.2). Furthermore, as discussed by Padovani et al. 
(2009) and Padovani (2011), R is useful for quasar samples, 
where it is related to the jet/disk ratio, but loses its meaning 
for radio galaxies. 


part of the LF and, when integrated over the full range 
of powers, the resulting /rlagn is much smaller. 

The best way to approach the question of the in¬ 
trinsic value of /rlagn is to look at the band under 
discussion, namely the radio one. We plot in Fig. 8 our 
best estimate(s) of the local LF for RQ (radiative-mode) 
AGN, together with the z < 0.3 LF of radiative-mode RL 
AGN from Best et al. (2014) (dotted line). One can see 
that the fraction of radiative-mode RL AGN is a strong 
function of radio power, as the shape of the two LFs is 
very different. In particular, for Pi. 4 GHz > 10 24 W Hz -1 
radiative-mode RL AGN dominates and /rlagn > 0.5, 
while for Pi. 4 GHz JS 10 22 W Hz -1 , which is the limit¬ 
ing value for radiative-mode RL AGN, /rlagn ~ 0.04 
(where we have used the double power law fit derived in 
the previous sub-section). To derive the overall intrinsic 
value one would need to integrate down to lower values 
but we do not yet have a handle on the low-end of the 
radiative-mode RL AGN LF. 

Our RQ AGN include absorbed (type II) and unab¬ 
sorbed (type I) sources, while the LF of radiative-mode 
RL AGN from Best et al. (2014) is limited to radio galax¬ 
ies and steep-spectrum radio sources (Sect. 6.1.1). How¬ 
ever, as radio quasars are relativistically beamed, they 
tend to be much more luminous and less numerous than 
radio galaxies and therefore we believe that down to 
Pi. 4 GHz ~ 10 22 W Hz -1 their exclusion makes very little 
difference. 

We note that if one takes as dividing line between RQ 
and radiative-mode RL AGN the power at which the two 
LFs cross, one gets Pi.4GHz « 10 24 W Hz -1 at z ~ 0, i.e., 
the same dividing value between the two classes that has 
been suggested in the past (see Sect. 2 of Padovani 1993). 
Given the strong redshift evolution of both LFs, one can 
predict that this power is redshift dependent and will get 
larger at higher redshifts, as indeed appears to be the case 
(Padovani 1993). Future radio studies of deeper/wider 
samples, which will include RL AGN of the jet-mode and 
radiative-mode type, should be able to study the RL frac¬ 
tion dependence on redshift by deriving the evolving LFs 
for both classes. 


7 SUMMARY AND CONCLUSIONS 

We have used a deep, complete radio sample of 680 ob¬ 
jects down to a 1.4 GHz flux density of 32.5 /rjy selected 
in the E-CDFS area to derive the number counts of the 
various sub-mjy classes and to study in detail the evolu¬ 
tion and luminosity functions of radio faint AGN up to 
z ~ 4. Our main results can be summarized as follows: 

(i) What we consider to be the best determination of 
the source population of the sub-mjy radio sky shows 
that star-forming galaxies and AGN make up a roughly 
equal part of the sub-mjy sky down to 32.5 /rJy, with 
the former becoming the dominant population only below 
~ 0.1 mjy. RQ AGN are confirmed to be an important 
class of sub-mjy sources, accounting for ~ 25% of the 
sample and ~ 60% of all AGN, and outnumbering RL 
AGN at <0.1 mjy (Sect. 3). 

(ii) The AGN that make up the faint radio sky consists 
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of two totally distinct populations: RQ AGN, mostly of 
the radiative-mode type, and RL AGN, largely of the 
jet-mode kind, characterised by very different evolutions, 
LFs, and Eddington ratios (Sects. 5 and 6). 

(iii) RQ AGN evolve in radio power as oc (1 + z) 2 ' 5 , 
similarly to star-forming galaxies, up to z ~ 4 but with a 
hint of a slowing down above z ~ 1.3. Their LF is steep 
(<F(P) oc P" 2 - 5 for P > 2 x 10 22 W Hz" 1 ) and their 
T/pEdd is typically > 0.01 (Table 3 and Sects. 6.1.2 and 
6 . 2 ). 

(iv) RL AGN evolve in number oc (1 + z) 2 but only up 
to z ~ 0.5, beyond which their number density declines 
steeply oc (1 + z) -4 . Their LF is flat (<1?(P) oc P” 1 ’ 4 ) and 
their P/pEdd is typically < 0.01 (Table 4 and Sects. 6.1.1 
and 6.2). 

(v) The first determination of the local radio LF of 
RQ AGN appears to be consistent with that of Seyferts. 
Putting the two together, one derives a number density 
of RQ AGN equivalent to that of X-ray selected AGN, 
which shows we are looking at the same population from 
two different bands (Sects. 5.3.1 and 6.1.2). 

(vi) We have approached the question of the intrin¬ 
sic value of the fraction of radiative-mode RL sources, 
/rlagn, within the AGN population by using, for the 
first time, the radio band (while all previous attempts 
were normally done in the optical band). By combining 
the local radio LF of Seyferts and our RQ AGN with 
a previously published estimate of the LF of radiative¬ 
mode RL AGN, we find that for Pi. 4 GHz > 10 24 W Hz -1 
RL AGN dominates and /rlagn > 0.5, while down to 
Pi.4GHz - 10 22 W Hz- 1 /rlagn - 4% (Sect. 6.3). 

(vii) The surface density of radio-selected, radio-quiet 
AGN, ~ 860 deg -2 , is already larger than that reached 
by the deepest optical surveys. This means that sub- 
mjy radio surveys, given the appropriate ancillary multi¬ 
wavelength data, have already the potential of detecting 
large numbers of radio-quiet AGN by-passing the prob¬ 
lems of obscuration, which plague the optical and soft 
X-ray bands (Sect. 6.1.2). 

In an upcoming paper (Padovani et al., in prepara¬ 
tion) we will discuss in detail the evolution and LF of 
the SFG in our sample and how they compare to those 
of RQ AGN, which is very relevant for the issue of the 
mechanism behind radio emission in RQ AGN (see also 
Padovani et al. 2011 and Bonzini et al. 2015, submitted). 
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Figure Al. Redshift versus .3.6/mi flux density for our sources 
with reasonable and secure redshift information (quality flag 
^ 2 as defined in Bonzini et al. 2012). The solid line is the 
best fit for z 0.3, while the dashed lines represent the root 
mean square error. See text for details. 
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APPENDIX A: REDSHIFT ESTIMATES 

As shown in Fig. Al, reasonable and secure redshifts 
(i.e., spectroscopic with quality flag > 2, as defined in 
Bonzini et al. 2012) are strongly correlated with f 3 .e 1 u.rn, 
albeit with some scatter in particular at low redshift. The 
best and simplest approach to estimate the redshift for 
the ~ 8% of the objects in the sample without observed 
redshift is then to derive it from their f 3 . 6 fim by using 
the relationship shown in the figure (solid line), that is 
logz = —0.677 logf 3 .eium + 1.185. This was derived ap¬ 
plying the ordinary least-square bisector method (Isobe 
et al. 1990), which treats the variables symmetrically, to 
the z > 0.3 sample. The low redshift cut was applied for 
three reasons: 1. in order to reduce the scatter in the in¬ 
put data; 2. because the maximum / 3 . 6 /jm of the sources 
without redshift is ~ 200 microjy, which corresponds to 
z ~ 0.3 deriving the best fit for the whole sample; 3. 
because the median f 3 .efj.rn for sources without redshift 
is ~ 10 microjy and therefore they should typically be 
at large (z > 1) redshift; inclusion of the low-redshift 
sources in the sample might therefore bias our estimates. 
We note that the best fits for the different sub-classes give 
redshift values, which differ only by ~ 0.4 for sources hav¬ 
ing f 3 . 6 fim ~ 10 microjy (the median for sources without 
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redshift). Note also that for / 3 . 6 /im = 1, the minimum 
value for the sample, 2 ~ 15.3, which is on the high side. 
We therefore set the maximum redshift equal to the max¬ 
imum observed values (spectroscopic or photometric) for 
the three sub-classes, namely z ~ 4.7, 7.0, and 4.5 for 
SFG, radio-quiet, and radio-loud AGN respectively. 


APPENDIX B: THE LARGE SCALE 
STRUCTURE OF THE E-CDFS 

Gilli et al. (2003) have studied the large-scale structure 
(LSS) in the CDFS in the X-ray band, finding two main, 
narrow (A 2 < 0.02) structures at 2 = 0.67 and 0.73 
and other spikes at 1.04, 1.62, and 2.57. Indeed, the red¬ 
shift distributions shown in Fig. 2 peak in the 0.5 — 0.75 
bin for all classes. Given the small area of our survey 
one could worry that such redshift spikes might influ¬ 
ence some of our results. We then studied the LSS of the 
E-CDFS by looking for narrow (Az ~ 0.01 — 0.03) struc¬ 
tures and smoothing the observed spectroscopic redshift 
distribution in large (Az ~ 0.1 — 0.3) bins to estimate 
the background distribution. Redshift peaks were defined 
as spikes being > 2<r above the background for Poisso- 
nian errors. Applying this procedure to the whole sample 
we recovered the main structures of Gilli et al. (2003), 
namely those at 2 ~ 0.67, 0.73, and 1.62. We also find a 
structure at 2 = 0.52, present in the near-IR band data 
of Gilli et al. (2003), and a new one at 2 ~ 0.12. Studying 
the LSS for the various sub-classes, however, showed that 
not all redshift peaks are present for all classes. We then 
studied the LSS individually finding the following signif¬ 
icant peaks: 1) SFG: 0.122 - 0.126 and 0.729 - 0.736; 
AGN: 0.521 - 0.537, 0.661 - 0.686, 0.720 - 0.738, and 
1.603-1.619; RQ AGN: 0.521-0.529, 0.652-0.671, and 
1.591 - 1.619; RL AGN: 0.675 - 0.686 and 0.731 - 0.734. 
We then created “no LSS” sub-samples by removing the 
sources in the bins listed above in excess of the back¬ 
ground distribution. The selection of the sources to be 
kept, which typically included ~ 20% of those in the bin, 
was done by picking the object (s) with radio flux densi¬ 
ties closest to the median value of the bin. 
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